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Abstract 

The manifestation of measurements, randomly distributed in time, on the evolution of quantum 
systems are analyzed in detail. The set of randomly distributed measurements (RDM) is modeled 
within the renewal theory, in which the distribution is characterized by the probability density 
function (PDF) W{t) of times t between successive events (measurements). The evolution of the 
quantum system affected by the RDM is shown to be described by the density matrix satisfying 
the stochastic Liouville equation. This equation is applied to the analysis of the RDM effect on the 
evolution of a two level systems for different types of RDM statistics, corresponding to different 
PDFs W{t). Obtained general results are illustrated as applied to the cases of the Poissonian 
[W(t) ~ e~ Wrt ] and anomalous [W(t) ~ l/t 1+a ,(a < 1)], RDM statistics. In particular, specific 
features of the quantum and inverse Zeno effects, resulting from the RDM, are thoroughly discussed. 

PACS numbers: 03.65.Xp 
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I. INTRODUCTION 



The effect of measurements on the evolution of quantum systems has recently been stud- 



ied very active 
example, refs. 



y_ both experimentally and theoretically (for comprehensive review see, for 
4]). The interest to this problem was inspired by the pioneering paper 
concerning the analysis of the specific feature of the manifestation of measurements which 
was called the Zeno effect, showing itself in the strong decrease of the decay rate of the 
quantum state with the increase of a number of measurements [5( • Since then a lot of works, 
analyzing different aspects of this effect, have been published. 

Traditionally the quantum Zeno effect is considered assuming a set of measurements 
to be distributed equidistantly in time (with constant time between measurements). The 
equidistant distribution results evidently in simplification of the mathematical treatment of 
the problem and experimental observation of the measurement effect IrM]. 

Naturally, the majority of manifestations of the measurement effect (including the Zeno 
effect), found for equidistant distribution of measurements, are expected to occur in the case 
of irregular distribution as well. However, some additional analysis is certainly needed. 

In this work we will discuss the measurement effect in the interesting special case of 
the irregular distribution, the case of measurements randomly distributed in time, called 
hereafter randomly distributed measurements (RDM). The random process of measurements 
is modeled within the renewal approach (RA) which treats the sequence of measurements 
as a stochastic set of renewals jf], Q]. In the RA the distributions of time intervals t between 
successive renewals are assumed to be stochastically independent and are described by the 
probability density function (PDF) W(t) (often denoted as ip(t) Q]). 

The evolution of the quantum system, affected by the RDM of this type, is shown to be 
described by the stochastic Liuvilles equation (SLE) for the density matrix of the system 
. The SLE allows one to analyze the effect of the RDM for different types PDFs W(t) 
in a fairly simple analytical form. 

In our work general expressions for density matrix of the system, affected by the RDM are 
derived. With the use of these expressions some important specific features the RDM effect 
are analyzed for different types of W(t) behavior. In particular, the quantum Zeno effect, 
i.e. the decrease of the decay rate of the state with decreasing the average time between 
measurements i, is predicted only in the case of rapidly decreasing PDFs W(t) for which the 
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average time between measurements i = J^drrW^r) is finite. As to heavy tailed PDFs, 
which cannot be described by a finite ?(?—>■ oo), for these PDFs the quantum Zeno effect is 
shown to be absent. Such a drastic difference of the RDM effect in these two cases appears 
to be conveniently treated in terms of the Zeno-effect dependence on the characteristic rate 
w r of the Laplace transform W(e) as a function of the Laplace variable e. The fact is that, 
unlike t, the rate w r can be introduced independently of the mathematical form of W(t) 
decrease, as it will be shown in our work, though for rapidly decreasing W(t) these two 
parameters are closely related: i ~ l/w r . 

To illustrate the general results we will discuss the case of the Poissonian distribution 
with W(t) = e~ Wrt , as an example of the RDM with rapidly decreasing PDF, and the case 
of anomalous heavy tailed distribution W(t) ~ l/t 1+a with a < 1. 

In particular, within the Poissonian model some characteristic properties of the RDM ef- 
fect are studied which are typical for case of rapidly decreasing PDF W(t). These properties 
manifest themselves in specific features of the time dependence of the probability p(t) to 
survive in the measured state. It is shown, for instance, that in the limit of small character- 
istic time t r = w^ 1 between measurements p(t) is the exponentially decreasing function with 
the rate non-monotonically depending on t r : the decreasing dependence at very small t r , 
corresponding to the quantum Zeno effect, is changed by the increasing one, associated with 
the inverse Zeno effect [2-4, 10|. as t r is increased. The time t Tm of the change is determined 
by the parameters of system. 

In both Poissonian and anomalous RDM models the manifestation of relaxation and, for 
example, the appearance and specific features of the quantum and inverse Zeno effects in 
the presence of relaxation are also analyzed be means the proposed approach. 

In conclusion, possible realistic examples of processes, in which the predicted effects can be 
observed, are thoroughly discussed. In particular, some realizations of both Poissonian and 
anomalous RDM in the processes with the participation of Brownian particles are considered. 

II. FORMULATION OF THE PROBLEM 

In this work we consider the effect of the RDM on the evolution of dynamic quantum 
systems, i.e. the systems, in which relaxation is absent. The statistical properties of the 

n q 

RDM are treated within the renewal theory |6|, LZj , well known in statistical physics. 
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In order to analyze the RDM effect we, first, should clarify some details of this effect for 
the finite number, say, n of measurements. 



A. Measurement affected quantum dynamics 

The dynamic quantum system is characterized by the wave function \if)(t) ) (the vector in 
the Hilbert space) whose time evolution is governed by the Schrddinger equation (h = 1) 

\ip) = -iH\ij)), with H = H Q + Hi. (2.1) 

In this equation H is the Hamiltonian of the system represented as a sum of the free and 
interaction parts, H and H iy respectively. The initial condition for eq. (12.1 j) \if)(t = 0)) = 
depends on the process considered (see below). 

The time evolution of the wave function is usually described by expansion in some (com- 
plete) basis of functions 1^) in the Hilbert space. In our further analysis, for the sake of 
convenience, we will assume to be basis of eigenfunctions of the free Hamiltonian Hq. 

In accordance with the conventional von Neumann rule 11], measurements, which 

show that the system is in some state suggested to be the eigenstate of H with the wave 
function \ip ) (H \ip Q ) = cu \i/; )), enable one to obtain the probability pit) to find the system 
in the state p(t) = |(^ |^W)| 2 - 

In what follows we will analyze theprobability p(t) of survival in the state \ipo) after a set 
of measurements. Following refs. jl, [J the subtle problem of the evaluation of the multiple 
measurement effect will be treated assuming that "if every time the measurement has a 
positive outcome and the system is found in the initial state, the wave function "collapses" 
and the evolution starts anew from l^o)"- I n this case the problem is to calculate the 
probability p(t) to find the system, initially created in (i.e. with \ipi) = \ipo)), i n the 
same state |^ ) after the measurements. 

It is worth noting that, in fact, the effect of measurements reduces to population and 
phase relaxation in the system under study (see below). In such a case, in general, the 
evolution of the system should be described with the density matrix p(t) satisfying the 
Schredinger like equation but in the Liouville space |ll|. In this space the matrix p(t) is 
represented as a vector in the basis, consisting of bilinear combinations of wave functions 
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\jj') = \^j)^A-- (2-2) 
p{t) = \p(t)) = J2 ..,M*)I">- (2-3) 

Note that in the Liouville space for any vector \v), {y — jj'), the corresponding conjugated 
one (u\ is defined by the relation (u\i/) = 8 w i. Noteworthy is also that vectors \jf) and 
\fj) are considered as independent and = 0. 

For the dynamic system, whose wave function \ip{t)) satisfies eq. (12. ip . the density 
matrix is represented as p(t) = \ip(t)){ip(t)\. The above mentioned Shrodinger equation (in 
the Liouvulle space) for this matrix is similar to eq. ( 12. ip but with the Hamiltonian operator 
H replaced by the superoperator H: 

p = -iHp, where Hp = [H, p] = Hp- pH. (2.4) 

According to the definition ( 12. 2 p the matrix elements of the superoperator H are given by 
(kk'\H\jj') = (ipk\H\il)j)8j'k' — (i>j'\H\iJ}k')fijk- It is important to note that for the Hermitian 
Hamiltonian H the superoperator H is also Hermitian (but in the Liouville space). 

The effect of a measurement on the evolution of the quantum system is conveniently 
determined in terms of the projection superoperator P m = |00)(00|, where 1 00) = \xjjo)(i[ , o\- 
In particular, after n successive measurements at times {tj}(j = 0, . . . ,n) ordered as t > 
t n > t n _i > ■ ■ ■ > ti > to — 0, the density matrix p(t) can be obtained with formula 

p(t) = UtMf>o> (n = 0,l,...), (2-5) 
in which t n = (t n , . . . , t Q ) is the vector of measurement times, Ut (t) = Uo(t — t ) and 

n 

KM = U (t - t n ) HiPMtj - (2.6) 

with Uo(t) = e~ tHt being the evolution superoperator (operator in the Liouville space), and 
/o = |-0o><-0o| = 100). 

With the density matrix the observable under study, i.e. the probability of survival in 
the state \ipo) after n measurements, is expressed as 

Pn (t) = (iPoWMo) = ^[PmUtJt)}, (2-7) 
where the trace is evaluated over the states \jf) in the Liouville space. 
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The probability p n {t) is known to depend not only on the parameters of the Hamiltonian 
H but also on the number n of measurements. The dependence is thoroughly analyzed in 
a large number of publications assum ing of fixed time interval r m = t/n between measure- 



ments (see, for example, reviews 21-14. 1 1 II] ) . Below we will extend the analysis assuming the 



multiple measurements to be a stochastic processes (in general non-Markovian) , which will 
be modeled within the RA js, 7]. 



B. Renewal approach 

1 . General formulas 

In the RA the times t n of events (measurements) are considered to be randomly dis- 
tributed in time with intervals r n = t n — t n _i between successive events described as inde- 
pendent random variables with the monotonically decreasing PDFs W(r n ), the same for all 
intervals. Note that W(t) is defined only for positive t, with W(t < 0) = 0. To completely 
characterize the statistics of events one also needs the probability P(t) that the interval 
between the successive events is greater than t: P(t) = f°°dt'W(t'), which is, naturally, 
normalized by the relation P(0) = 1. 

In what follows we will mainly operate with the Laplace transforms denoted as 

POO 

Z(c) = / dtZ(t)e~ et (2.8) 
Jo 

for any function Z(t). In particular, noteworthy is the relation -P(e) = [1 — W(e)]/e and 
suitable representations 

W(e) = [1 + ^(e)]- 1 and P(e) = [e + e/$(e)]- 1 . (2.9) 

in terms of the auxiliary function $(e) 

The analytical form of $(e) is completely determined by that of W(t). In what follows, 
to specify the characteristic scale of $(e)-dependence we will introduce the characteristic 
rate w r whose meaning will become clear from some particular examples of functions $(e) 
considered below. It is important that this rate can be introduced both for rapidly and 
anomalously slowly decreasing PDFs W(t) (see below). 
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The effects analyzed in this work are essentially controlled by the small e behavior of 
$(e). In general, we can only say that $(e) e — >° 0. Some additional information on $(e) 
behavior at e — > can be obtained by the analysis of the long time dependence of W(t): 

a) If at t — > oo the PDF W(t) decreases so rapidly that average time (first moment) 
t = J^°dttW(t) ~ w~ x is finite, at small e the function $(e) can be represented as 

$(e) e =° ei+o(e/w r ). (2.10) 

b) If, however, W{t) is a heavy tailed function: W(t — > oo) — > l/t 1+a with a < 1, and, 
therefore, t is infinite (does not exist), then $(e) ~ (e/w r ) a jy, 7]. 

Within the RA, due to independence of renewals, the PDF Wt n (t) of n events (n > 0) 
at times {tj}(j = 0,1, ... ,n), satisfying the relation t n > t n -\ > ■ • • > t\ > to = and 
combined into the vector t n — (t n , . . . , ti), is given by 

n 

These functions completely describe the stochastic renewal process. In particular, with 
the use of W tn the probability ir n (t) to observe n events in the time interval (0, t) is expresses 
by fl: 

7T n (t)= / df P(t-t')W n (t')dt', (n>0). (2.12) 
Jo 

In this formula W n (t) is the PDF of n events, which for n > 2 is equal to Wt n integrated 
over all tj<„_i (from to tj + i, respectively) except for t n = t, and for n = 0, 1 are defined 
by W\{t) = W(t), and Wo(t) = 5(t). The functions n n (t) can, evidently, be represented as 

WJt) = / deW n (e)e et de. (2.13) 

2vrz J_ lOD 

Note that according to eq. ( 12.1 3p W (t) = 5(t) and therefore 7r (t) = P(t). It addition, it is 
worth noting, that the probabilities n n {t) satisfy the normalization condition Xl^lo 71 ™ = ^■ 

2. Examples of renewal processes 

In our further analysis special attention will be paid to some particular renewal processes 
corresponding to different distribution of times of renewals, i.e. different functional form of 
the PDF W(t). 
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a. Poissonian distribution. The most simple is the Poissonian model, corresponding 



to 



BE 



P(t) = e~ Wrt , W{t) = w r e~ Wrt , (2.14) 
so that W(e) = [1 + (e/uv)] -1 and therefore 

$(e) = e/w r . (2.15) 

b. Equidistant distribution. The model of equidistant distribution, in which 

P(t) = d(t r -t), W(t) = 8(t-t r ), (2.16) 

with t r = l/w r , describes the set of events with the constant interval t r between successive 
ones [12]. This is just the model considered in almost all publications concerning the Zeno 
effect. In this model 

$(e) = e tre - 1. (2.17) 

c. Anomalous distribution. The anomalous model implies the Levy-type distribution 
[l3| of times between events with the heavy tailed behavior of W(t) ~ l/t 1+a , (a < 1). 
One can find a number of such type of models which predict the same results for long time 
features of processes. In our analysis we use the simple one, for which 6 s| 

P(t) = E a (-(w r t) a ), W(t) = -P(t), (a < 1), (2.18) 



where E a (-x a ) = (2m)- 1 fj^dz e xz (z + z 1 ^)- 1 is the Mittag-Leffler function [13|. This 
model corresponds to 

$(e) = (e/w r ) a . (2.19) 
It is easily seen that the expression (12. 18)) predicts the heavy tailed behavior the PDF W(t). 



III. STOCHASTIC LIOUVILLE EQUATION 

The above consideration shows that the RDM effect on the evolution of the quantum 
system under study is expressed in terms of the the superoperator U(t) [eq. ( 12.61) ] averaged 
over the stochastic process of measurements. In what follows this averaged superoperator 
will be denoted as U(t). 
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The problem of finding U(t), which is a functional of type of ( 12. 6j) averaged over the 
renewal process, has already been discussed in literature [8|, |9j. It is shown to reduce to 
solving some equation for this operator, which is called the SLE. 

For some types of renewal process the averaging is essentially simplified. The most well 
known example is the Poissonian process, simplification for which results from the Markovian 
nature of this process. In the RDM analysis the Poissonian process of quantum transitions 
(jumps) is just a realization of the Markovian process of "migrations" (or jumps) over states 
\jj') of the quantum system under study (in the Liovulle space). The treatment of this 
Poissonian "migration" process within the widely accepted continuum time random walk 
approach [l^. [l5[] shows that the PDF pu for the "migrating" system satisfies the Markovian 
equation p M = —w r (l — P m )PM- The effect of such "migrations" on quantum evolution is 
described by the SLE for the evolution suproperator U (t) = —iHU (t) — w r (l — P m )U(t) [l^ 
(see Sec. IV.C). 

Fortunately, the SLE can be obtained for any type of renewal process. In the most_general 
form the SLE is rigorously derived with the Markovian representation of the RA jg, |9j. In 
this work, however, we will restrict ourselves to the simplest variant of the RA, for which 
the SLE can be obtained fairly easily. Below we will outline some details of the derivation. 

The derivation is based on the fact that in the absence of measurements (n = 0) the 
evolution operator Ut (t) = Uo(t), while for any number n > 1 of measurements the operator 
Ut n (t) [eq. fl2.6p ] and the PDF of measurements W tn (12.111) are represented as products of 
terms depending on differences of times tj—tj-i. In such a case we get the following formulas 
for corresponding average evolution operators U n (t): Uo(t) = P(t)Uo(t), for n — 0, and 

V n (t) = [ dt n P(t-t n )W u U u (t), forn>l, (3.1) 
Jo 

where dt n = YYj=idtj- The convolution-like form of formulas for U n (t) results in a simple 
representation for the Laplace transforms of these functions: 

U n (e) = P(Cl £ ) [P m W(Cl e )] n with (l € = e + iH . (3.2) 

Summing up the contributions for different numbers of measurements we thus get the 
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evolution operator U(t) averaged over the renewal process |8j, |9] 

~ 00 ^ 

U(e) = U n {e) = P(Cle) [1 - P m W{Cl e )] 1 (3.3) 

n=0 

= ^(fiJ^^ + Qj]" 1 . (3.4) 

where Q m = 1 — P m = J2jj>^oo \ jf) Uf \ * s the superoperator of projection onto the subspace 
{././)}• (././' / ('()). 

The expression (13. 4p can be treated as a solution of the non-Markovian SLE for U(t) 

JJ = -iMJ -Q m [ dr M(T)e~ iilT fj{t - r) (3.5) 
Jo 

in which 

1 c\ r°° 

M(t) = — - dee^-\e). (3.6) 
2m dt J_ ioo 

is the memory function whose analytical properties are essentially determined by those of 
the PDF W(t). 

The SLE (13.51) is not quite convenient for applications. In our further analysis we will 
mainly use the expression ( 13. 4p for the Laplace transform of the evolution operator. 

In accordance with the general formula (12.71) the (average) probability p(t) of survival in 
the initial state (l^o)) is completely determined by the average evolution operator U(t). In 
terms of the Laplace transforms the corresponding expression is written as 

p(e) = Tr[pj7(e)] = Tr[PjJ(e)P m ], (3.7) 

where, similar to formula (12. 7p . trace is evaluated over the states \jj') in the Liouville space. 

In the conclusion of this general analysis we will show that with the use of above- 
obtained formulas p(e) can be expressed in terms of the probability po(t) = (00\Uo(t)\00) = 
(00|e~ im |00) of survival in the initial state |^o) m the absence of measurements. According 
to eq. (13. 7p the Laplace transform p(e) is determined by the trace of the supermatrix 



U P (e)=P m U(e)P m = J2 UpM (3.8) 
where U p n (e) = P m U n (e)P m . Each term U p„{e) in the sum can be represented by formula 

F P Je)=Pp(e)%e). (3.9) 
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In this expression we have introduced the supermatices 

X P (e) = P m X{(l e )P m for X = P,W, (3.10) 

which are directly related to the probability po(t), as is clear from the relations 
P m X(tt € )P m = J °° dte- et P m e- iflt P m X(t) and P m e~ i{lt P m = P m p {t). Finally we obtain 

U Pri (e) = P m P Po (e)W p n (e), (3.11) 

where 

POO 

X P0 (e)= dte- et X(t)p (t) for X = P, W. (3.12) 
Jo 

Substitution of formula ( 13. lip into eqs. ( 13. 8ft and ( 13.71) yields 

p(6) = P P0 (e)[l-%,(e)]- 1 . (3.13) 

In principle, both formulas (13.71) [with eq. (I3.4j) ] and (13 . 13[) are quite suitable for the 
analysis of the RDM effect on p(t). In what follows, however, we will mainly apply the 
formulation based on the first formula [( 13.71) ]. though, eq. ( 13.131) will also be used to clarify 
some particular properties of the effect. 

IV. ZENO EFFECT ON TWO LEVEL SYSTEMS 

In this section we will analyze the effect of the RDM on some model quantum system to 
illustrate the specific features of manifestation and treatment of the Zeno effect in this kind 
of measurements. 

The simplest (though quite realistic) system, which enables one to significantly simplify 
mathematical problems, is the two level system. It is very convenient for detailed description 



of all important features of the Zeno effect 



A. Hamiltonian of the two level system 

In our analysis we will use the Hamiltonian in the form [see eq. ( 12. ip ] 

Ho = - |2)(2|), H t = v(\l}(2\ + |2)(1|). (4.1) 

in which e and v are a positive real parameters. The measured state is assumed to be the 
state |1), i.e. |^ ) = |1). 

11 



In the Liouville space the Hamiltonian is represented as a 4 x 4 matrix 



H 



2e(|12)<12|-|21)<21|), 



(4.2) 



Hi 



v[(\U) - |22))«21| - (12|) 



+ (|21)-|12))((11|-(22|)], 



(4.3) 



in the basis \jf) = [defined in eq. fl2.2jl ] . In this basis the superoperator P m , describing 



the effect of a measurement, is written as 



P m = |11)(11| and Q m = V \jf)(jfl (4.4) 



The Laplace transform of the function under study, the survival probability p(e), can 
conveniently be evaluated with the expression (13.71) . which, as applied to the two level 
system considered, is written as 



With the above formulas at hand one can analyze the specific features of the Zeno effect 
for any type of the RDM. 

B. General results 

Here we will obtain some general results, valid for quantum systems with arbitrary (but 
finite) number of levels, to clarify the manifestation of the analytical properties of the de- 
creasing function W(t) in the specific features of the Zeno effect. 

The most important property of function W(t) is the rate of decrease at long times which 
is determined by the behavior of $(e) in the limit e — > (Sec. IIB.l), i.e at e w r , where w r 
is the above defined rate, characterizing the time t r — l/w r of the onset of the asymptotic 
long time behavior of W(t). Just in this limit, or more accurately in the limit of large 
w r , when £ = ||f2 e ||/w r <C 1, one can demonstrate some important properties of the RDM 
effect. In principle, essential conclusions can be made without any particular assumptions 
on $(e)-behavior except for the relation $(e) e — *-° 0. However, in the analysis it is suitable 
to keep in mind the approximation $(e — > 0) ~ (e/w r ) a , with a < 1, which is of particular 
interest for our further discussion. 



p(e)=Tr[P m f/(e)] = (ll|f/(e)|ll>. 



(4.5) 
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The above mentioned important general conclusions concern the properties of p(e) at 
£ <C 1. In this limit we get from eqs. H3.4[) and (14. 5 p 

p(e) « Poo(e) = Tr[P m fi7 1 $(a)]/(ll|$(a)|H) 

= (lii^^ii^/aii^^in). (4.6) 

This general expression is valid for systems with any number of levels, though it is presented 
in terms of the considered two level model (I4.2p .( l4~3j) . 

Formula (I4.6P can be derived from eqs. (13. 4p and (14. 5p by taking into account some 
characteristic properties of the matrix L = $(fi e ) + Q m - For brevity, in our further study we 
will use the notation <&(f2 € ) = $. To clarify the derivation we will, first, analyze the specific 
features of two parts of L: L P = P m LP m = P m &P m and Lq = Q m LQ m = Q m &Q m + 
Q m , operating in subspaces {|11)} and {\jf ^ H)}? respectively. The fact is that in 
the considered limit £ = ||f2 e ||/w r <C 1, corresponding to ||$|| <C 1, the eigenvalue Xp = 
(11|$|11) <C 1 of the matrix Lp is much smaller than Lg-eigenvalues Xq v ~ 1, whose 
magnitudes are mainly determined by Q m . These estimations yield for the characteristic 
splitting of eigenstates 5X = Xq v — Xp the value SX ~ 1. 

The obtained splitting appears to be much larger than the L-induced interaction Lqp = 
PmLQm + QmLP m = P m $Q m + Q m $P m between the states of {|11)} and {\jf ^ 11)} 
subspaces: ||-£qp|| ~ ||$|| SX. This means that in the leading order in ||$||/5A ~ ||$|| <C 1 
the eigenstates and eigenvalues of the matrix L coincide with those of Lp + Lq. Of special 
importance is the coincidence of the lowest eigenstate \l) of L with |11): |/) ~ |11) + 
(jj'\jf) with Qji ~ ||$|| 1, which results in formula {l\L\l) = (ll|Lp|ll)[l + 
0(||$||)] with (ll|Lp|ll) = (11|$|11). This formula implies that with the accuracy ~ 
|| $11 <C 1 the Green's function L _1 is mainly determined by the contribution of the eigenstate 
10 -ill): 

L- 1 = P m (ll|<f>|ll)- 1 [l + 0(||l>||)]. (4.7) 

Other eigenstates (which with high accuracy coincide with those of Lq) make a contribution 
much smaller than that of |/), since the corresponding eigenvalues Xq v ~ 1 are much larger 
than X P = (11|$|11) < 1. 

Substitution of the expression (14.71) into eqs. (13. 4p and (14. 5 p yields formula (14. 6p . Note 
that in the same way it can also be derived by means of eq. (I3.13p . 
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This formula allows one to analyze the limiting behavior of the RDM effect in the limit 
w r — > oo, i.e. for very small characteristic time between measurements. 

In particular, we can discuss the case of rapidly decreasing W(t) for which $(e) ~ ei 
[see eq. f l2.10p ]. Taking into account the relation (I2.10p one arrives at the estimation 
p(e) = e _1 [l + 0(||f2||/uv)] which means that 

p(t) w i^° i. (4.8) 

The limiting relation (14. 8[) demonstrates the localization of the system in the measured 
state in the limit w r — > oo, which is associated with the quantum Zeno effect (2-4]) expressed 
in terms of the RA. 

The dependence of p{t) on the average time i ~ w~ l between measurements is, in general, 
not universal. It is also determined by other parameters of the system, for example, the 
parameters of the Hamiltonian. In such a case, to characterize the quantum Zeno effect as 
a function of w~ l , one can apply the average time tz{w r ) = f™dtp(t) = p(e = 0): 



e-i-0 



= P Pi m-W Pl (e)r\^ - (4-9) 

where Q e = e + iH. The second of eqs (14.91) is written with the use the expression (13. 13ft 
for p(e = 0) in which, however, the functions P Po (e) and W Po (e) are replaced with P pi (e) 
and W Pl (e), respectively, defined by X Pl (e) = j^dt e" e 'X(t)pi(t) for X — P,W. This 
replacement is made due to the change of notation for the measured state (|1) instead of 
|0)), according to which, to avoid possible confusions, po(t) should be replaced by p\{t) = 
(ll|e-^|ll). 

In general, tz{w r ) can be calculated only numerically. Some conclusions on the specific 
features of this parameter, however, can be made by the analysis of limiting behavior of 
t z {w r ) at w r — > 00 and w r —> 0. 

It is seen from the above definition of P Pl (e) and W^e) that the value of tz is finite 
in the case of finite i = J °° dttW(t) = J °° dtP(t). In general, tz can be calculated only 
numerically. However, some conclusions on the specific features of this parameter can be 
made by the analysis of limiting behavior of tz(w r ) at w r — > 00 and w r — > 0. 

1) For w r — > 00 we get the relation tz{w r — > 00) — > 00 which follows from the definition 
of this parameter and the relation p(e) = e _1 valid in this limit. 
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2) In the opposite limit w r — > the parameter tz{w r ) also grows to infinity: tz(w r — > 0) ~ 
l/w r oo. This dependence can be obtained by analyzing ^-dependence of P Pl (e = 0) 
and W pi (e = 0) at small w r . The fact is that the probability Pi{t) = (11 |e — 1 11) , which 
determines the values of these two functions, can, in general, represented as: Pi{t) = pi + 
5pi(t), where pi is independent of time and 5pi(t) is the oscillating part represented as a sum 
of harmonically oscillating functions. The most important for our analysis is pi, which can be 
found by expansion of the evolution operator e~ lHt in the basis of eigenfunctions \<pj) of the 
Hamiltonian H: pi = KV'ilv^i)! 4 < 1- With this representation for pi(t) one can obtain 
the estimations P pi (0) ~ pi J^°dtP(t) = pit = pi/w r and W P1 (0) ~ pi J^°dtW(t) = pi < 1, 
in which the contribution of the oscillating part Spi(t), negligibly small in the limit w r — > 0, is 
ignored. Substitution of these relations into eq. (14. 9 p leads to the above-mentioned limiting 
dependence tz(w r — > 0) ~ pi(l — pi)^ 1 w~ 1 . 

The above analysis shows that the behavior of tz{w r ) is non-monotonic with the minimum 
of this function at some w r whose value is determined by the parameters of system. The 
validity of such a conclusion will be demonstrated below as applied to the case of the 
Poissonian RDM distribution, as an example. 

C. Poissonian distribution of measurements 

In the case of the Poissonian RDM statistics, when $(f2 e ) = Cl e /w r , equation (13.41) for 
;he Laplace transform U(e) is essentially simplified reducing to the SLE of Schrodinger type 



which is briefly discussed in Sec. III. 

For the two level system this equation is the system of four linear equations which can 
be solved analytically. We are not going to present the cumbersome expressions for matrix 
elements of U(e) but restrict ourselves to obtaining only one element corresponding to the 
observable under study, the survival probability p(t) [see eq. (j2.7j) ] . For the initial condition 
l^o) = |1) the Laplace transform p(e) is given by 




(tt e + w r Q m )U(e) = (e + iH + w r Q m )U(e) = 1. 



(4.10) 



p(e) = (H|t/(6)ill) 

e + w r + w(e) 



(4.11) 



e 2 + e[w r + 2w(e)] + w r w(e) 
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where 

w(e) = 2(e + w r )v 2 /[(e + w r ) 2 + 4 £ 2 ]. (4.12) 

Of special interest is, naturally, the limiting variant of general formula (14. lip correspond- 
ing to large average rate w r of repetition of measurements. It is seen from this expression, 
that in the limit w r ^> v (and for e < w r ) p(e) ~ (e + Wo)" 1 , i- e - for t > w' 1 , 

pit) « e"* * with w = w(0) = -f— . (4.13) 

+ 4e z 

Note that in the considered limit uv 3> u the rate wq ~ uv(t> /w r ) 2 <!« r 

Equation (14. 12j) shows that in the limit of large w r the decay of the survival probability 
turns out to be exponential with the rate Wq,, in principle, non-monotonically depending on 
w r . The non-monotonic behavior, however, can can correctly be described by eq. (I4.13P 
only in the case e ^> v, when within the region of validity of this formula (w r ^> v ) there 
exists the subregion of w r values, e 3> w r ^> v, in which v)o(w r ) is the increasing function of 
w r : wo(w r ) ~ w r . The non-monotonic dependence Wo(w r ) can be treated as the acceleration 
of p(£)-decay by measurements at relatively low measurement rates w r <C e (associated 
with the inverse Zeno effect) followed by the deceleration at large w r ^> e: Wo(w r ) ~ l/w r 
(corresponding to the quantum Zeno effect). 

In general, the specific features of the RDM effect can be demonstrated with the param- 
eter tz(w r ) = p(0). The expression (14.111) yields simple formula for this function: 

t z = w; 1 + Wo 1 = iT^l + \{w 2 r + 4e 2 )}/w r (4.14) 

with e = e/v and w r = w r /v. In agreement with general qualitative conclusions (Sec. IVB), 
tz(w r ) non-monotonically depends on w r with tz ~ l/w r —too,tz ~ w r /{2v ) — > oo, 
and the minimum at 

w rm = (2 + 4e 2 y/ 2 . (4.15) 

For e > t; in the limit w r 3> v the non-monotonic behavior of tz(w r ) ~ w^ 1 with w Tm = 2e 
indicates the occurrence of two regimes of the Zeno effect, quantum and inverse, mentioned 
above. It is seen that in the case e ^> v the coordinate w rm {= 2e) predicted by eq. (14. 15ft 
coincides with that of the maximum of the rate v)o(w r ). The results of this analysis demon- 
strate that the parameter tz(w r ) is certainly useful for studying qualitative specific features 
of the RDM effect. 
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Characteristic properties of the behavior of the probability p(t) = p{r\r r ) (at fixed r = tv) 
as a function of the dimensionless average time r r = w' 1 = v/w r between measurements are 
shown in Fig. 1. 

At very small r r the probability p{r\r r ) monotonically approaches 1 as r r — > 0, in agree- 
ment with predictions of the quantum Zeno effect. The above analysis with the SLE (I4.10p 
clearly reveals the mechanism of the slowing down of 1 — > 2 transitions by the RDM. Accord- 
ing to this equation the RDM effect on the evolution of the state |1) is actually equivalent to 
the effect of the decay of the state |2) with the rate w r , which is accompanied by dephasing 
(decay of the density matrix elements (l|p|2) and (2|p|l)) with the same rate w r . Just fast 
dephasing in the limit w r — > oo leads to the strong reduction of the 1 — > 2 transition rate 
[~ v (v/w r )], associated with the quantum Zeno effect. 

At intermediate values of r r the function p(t\t t ) non-monotonically depends on r r with 
the minimum at some r rpm = v/w rpm . The position of the minimum is reasonably accurately 
estimated by means of w Tm (14. 15j) for all values of e = e/v used in numerical calculation: 
r rpm = vr rpm « w~\ The accuracy of this estimation is especially good for largest considered 
value e = 2.5, as expected from eq. (I4.13P and above analysis. 

The analysis of the expression (I4.13P makes it possible to understand the reason of the 
appearance of this minimum. It results from interplay between quantum oscillations in 
the system and RDM induced dephasing. Just this interplay manifests itself in the factor 
w 2 + Ae 2 in the denominator in formula (14.131) for wq, which is responsible for the non- 
monotonic behavior of p{r r ). 

In the region of validity of eq. (14. 131) . i.e. for r r = v/w r 1, the kinetics of p{r) decay is 
exponential. This means that the non-monotonic behavior of the dependence p(r r ) at a fixed 
t = vt stems from that of the rate uio(w r ). Such a dependence of the rate is traditionally 
treated as a transition from the quantum Zeno effect [in the region of decreasing behavior 
of wo(w r )} to the inverse Zeno effect [when wo(w r ) is the increasing function] [3]. In general, 
however, for non-exponential p(r)-dependence, it is not quite correct to associated the non- 
monotonic behavior of p(r r ) with any specific phenomena because of a large number of 
possible peculiarities of this behavior p(r)-kinetics in general. 



17 



D. Anomalous distribution of measurements 



Of special interest is the case of the anomalous RDM, in which the PDF W(t) anomalously 
slowly decreases at large times t: W(t) t ~^° l/t 1+a with a < 1. It is easily seen that this PDF 
cannot be characterized by the average time i between measurements in its usual definition 
though, of course, this functions still has the characteristic decay time (see below). In our 
further consideration we will apply the simplest Mittag-Leffier model for the anomalous 
PDF W(t), defined in eqs. (12.181) and (I2.19p . The existence and qualitative definition of the 
characteristic time t r = w' 1 in this model is clear from eq. (I2.19p . 

For brevity in our analysis of the anomalous case we will restrict ourselves to the most 
interesting limit of high characteristic inverse time between measurements w r , corresponding 
to ||fi||/w r <C 1. In this limit p(e) ~ Poo(e) can be evaluated with formula (I4.6p . which in 



the applied Mittag-Leffier model is represented as [17J] 



p^e) = (ll|fir 1 |ll)/(ll|^|ll) 
(2e 2 + 1) e"" 1 + fi?- 1 
~ (2e 2 + l)e a + n? ' 

where e — e/v and 



(4.16) 



QP = \ [( e + liv^WTiy + (e - 2ivV^TlY] . (4.17) 

It is easily seen that for the Poissonian RDM (a = 1) Poo(e) = so that Poo(t) = 1 as 
predicted in presence of the quantum Zeno effect. 

As for the anomalous RDM, the existence of the non-trivial limit Poo(e) itself indicates 
the violation of the charactristic Zeno-effect behavior of Poo( e ) i n this case. In the limit 
||fi||/uv 1 the function Poo( e ) [ an d thus p(t) = Poo(t)] appears to be independent of w r 
and determined by the parameters of the Hamiltonian only. 

At short times t e^.v^ 1 formula (14.161) predicts no transitions, i.e. Poo(e) ~ 1/e and 
pit) ~ 1. In the limit of long times t ^> e -1 ,u -1 , however, one gets 

Poo(t) ~ A a (t)r a } (4.18) 
where A a (t) is the oscillating function of time: A a (t) = a a + c a cos(2Et) + s a sin(2Et), 



in which E = v \/ e 2 + 1 and a a , c a , and s a are the constants depending on a and e. This 
asymptotic expression can be derived by taking into account that the most slowly decreasing 
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(and additive) contributions to the integral of the inverse Laplace transformation Poo(t) = 
(27ri)~ 1 Ja rfepoo(e)e e ' come from singularities of Poo(e) (14.16j) located at the imaginary axis 
of the complex plane of e. The singularities (brunching points) are determined by the terms 
e a ~ l and f2° _1 in the numerator in eq. (14.16)) . In the long time limit t ^> \j\Je 2 + v 2 these 
singularities contribute independently and the evaluation of contributions reduces to the 
calculation of integrals of type of {2ni)~ l dee a ~ 1 e et ~ t~ a , which leads to eq. ( 14. 18)) . 

It is important to note that the transition from the anomalous case a < 1 to the Pois- 
sonian one a = 1 is fairly non-trivial. In order to clarify the details of this transition we 
will consider the case of a close to 1, when 5 a — 1 — a -C 1 In this limit the onset of the 
inverse power type kinetics (14. 18ft displaces to very long times and the major part of the 
kinetics p(t) reduces to the exponential one, which can be obtained by expanding Poo(e) in 
powers of 5 a . For example, taking into account that for 5 a <C 1 x a ~ l w x(l — 8 a \nx), one 
can represent the numerator as (2e 2 + l)e a_1 + H" -1 = 2{e 2 + 1)[1 + 0(6 a )]. As for the 
denominator, similar expansion results in {2e 2 + 1) e a + f2° = 2(e 2 + l)(e + w z )[l + 0(5 a )}. 
Finally one gets Poo(e) ~ (e + Wz)' 1 and therefore 

Poo(t) w e~ Wzt with w z = |tt(1 - a)vVe 2 + 1. (4.19) 

It is worth noting that the characteristic behavior (14.191) is determined by that of Poo( e ) & t 
small |e| ~ 5 a . In this region the terms of higher order in 5 a = 1 — a, whose summarized 
contribution is denoted as 0(5 a ), result in the correction of the kinetics ~ 5^ \a(l/5 a ) which 
can be neglected in the limit 8 a <C 1. These terms, non-analytical in e, are responsible for 
the inverse time behavior of Poo{t) at long times. 

Figure 2 shows the dependence Poo{t) for some particular sets of parameters of the model. 



V. EFFECT OF RELAXATION 



A. General remarks 



So far we have considered the RDM effect on dynamical systems only, though it is known 
that the relaxation can strongly modify the manifestation of this effect |2|-|4|). In particular, 
in addition to the quantum and inverse Zeno regimes of the effect, analyzed above, there ap- 
pears another one, in which no strong influence of measurements on the evolution kinetics is 
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observed in the limit of short average time t ~ w r 1 — > 0. This regime can also be considered 



manifestation the inverse Zeno effect 



4| though this is a matter of definition. 



The influence of relaxation on the Zeno effect has already been analyzed within the model 



of equidistant ly distributed measurements 2^ J]). Here we will discuss this influence in the 
case of the RDM. 

In general, the analysis of the relaxation in quantum systems is a difficult problem. In this 
work we will restrict ourselves to the simple Markovian case, in which the kinetic equation 
for the density matrix can be written as 

p = — Cp, where C = iH + R. (5.1) 

Here R is the relaxation superoperator which can be evaluated in the short correlation time 
(or the Bloch-Redfield) approximation |18| . Some example of the relaxation model, leading 
to the particular expression for the superoperator R will be considered below. 

It is easily seen that in the presence of relaxation all general formulas are similar to those, 
obtained above for dynamic systems. The only difference is in the definition of Q e . In the 
presence of relaxation 

tt e = e + £ = e + R + iH. (5.2) 

To apply the results obtained in Sec. Ill with the redefined Q e we need to specify the 
relaxation superoperator, i.e. the relaxation mechanism. 

B. Simple relaxation model 

In order to illustrate the main features of the manifestation of relaxation in the RDM 
effect it is sufficient to analyze the simplest relaxation models. In our analysis we will discuss 
the variant of the model, widely accepted in the magnetic resonance theory in which 

R = w d (|ll)-|22»((ll|-(22|) 

+w p (|12)(12| + |21)(21|). (5.3) 

This operator describes population relaxation with the rate Wd and dephasing with the rate 
w p . The rates satisfy the relation w p > which ensures positivity of the density matrix 
pit) during the evolution, described by eq. (15. IH fl9 |. 

The model (I5.3P allows one to analyze fairly easily the important properties of the relax- 
ation effect for any type of the PDF W(t). 
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1. Poissonian distribution of measurements 



As we have already mentioned above, the case of the Poissonian RDM can be considered as 
a simple example of renewal processes with rapidly decreasing PDF W(t). The substitution 
of Cl t ( 15. 2 p into eq. ( 14. lUp leads to the expression (14.1 ip for p(e), but with w(e) replaced by 
w d (e) = w(e) + w d and with w r changed by w r + w p in the function w(e) itself: 

= e + w r + w d (e) 

e 2 + e[w r + 2w d {e)\ + w r w d (e) 

e + w d (e)(e + w r )/[e + w r + w d (e)Y ( 

where 

w)d(e) = w d + 2(e + w w )v 2 /[(e + Wrp ) 2 + Ae 2 }. (5.5) 

and Wrp = w r + w p . 

This expression enables us to analyze all possible effects of relaxation in the case of 
Poissonian distribution. Formula (15.41) predicts different types of p(e) dependence [and thus 
p(t)] on the average time between measurements t = w' 1 . The important parameter, which 
essentially determines the dependences of p(e) and p{t) on w r , is 

w° d = w d (0) = w d + 2w rp v 2 /(w 2 rp + As 2 ). (5.6) 

In particular, the parameter tz is directly related to w d : 

t z = p(0) = l/w r + l/w° d . (5.7) 

Here we summarize some most interesting limiting p(t) dependences on w r . 

a. Slow population relaxation, w d -C w r . In the case of slow population relaxation (or 
fast repetition of measurements as it is considered in Sec. IV C), when w d , w r (implying 
that v, w d w r ), formula (15. 4p strongly simplifies predicting exponentially decreasing pit): 

p(e) w (e + w^y 1 and p[t) w e~<\ (5.8) 

for which tz = ^-/w d . Depending on the relation between parameters of the system these 
expressions predict different m r -depedences of the RDM effect: 

1) For fast dephasing, when w p ^> w r and w rp ~ w p , we find that the p(t) decay rate and 
tz are independent of w r] j.e. of the measurements, which can be considered as the onset of 
the inverse Zeno regime 



w d « w d + 2v 2 w p /(w 2 + As 2 ). 
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2) For slow dephasing, when w p <C w r , the decay rate w° ~ Wd + 2v 2 w r /iw 2 + 4e 2 ) 
and therefore both quantum and inverse Zeno regimes are possible. For w r > e we get 
the decreasing function w^(w r ) ~ Wd + 2v 2 /w r , corresponding quantum Zeno case, while in 
the opposite limit w r < e we obtain the dependence w^(w r ) ~ Wd + (v 2 w r )/ (2e 2 ), which is 
associated with the inverse Zeno effect. Note that the above analysis of w®(w r ) dependence 
(and RDM effect, in general) is closely related to that of the rate Wo(w r ) of the (exponential) 
p{t) decay in the absence of relaxation and for fast repeated measurements (w r ^> v), which 
is described by eq. (14.131) . 

b. Fast population and phase relaxation Wd, w p ^> w r . In the limit of large population 
and phase relaxation rates Wd, w p ^> w r (recall that w p > Wd/2) of special interest is the case 
w° > w r , in which the evolution kinetics consists of two stages: the stage of fast equilibration 
(at t ~ and the stage of slow quasiequlibrium evolution affected by measurement (at 

t > l/w r > 1/w®). During the first fast stage the survival probability pit) decreases from 
1 to 1/2 according to pit) ps |(1 + e~ 2w ° dt ). After that (during the second most interesting 
stage) pit) decreases exponentially: 

pie) w |(e + Wr/2)- 1 and p(t) w \e~ {w ^ 2) \ (5.9) 

so that tz = l/w r in agreement with the prediction of eq. (15.7ft . Formulas (15. 9 ft can also 
be derived with the use of eq. (13. 13[) in which P Po (e) and W Po (e) should be replaced by 
P pi (e) and W pi (e), respectively, as it was mentioned in the analysis of eq. (14. 9p . Taking 
into consideration that in the limit of fast relaxation after short time ~ w® the survival 
probability Pi(t) = (l|p(t)|l) in the absence of measurements is given by pi(t) ~ 1/2, 
we get P Pl (e) = j^dt e - rf P(t) Pl (t) w (e + w r )- 1 and W Pl ie) = \ f™dt e~ tl W i^t) w 
\w r ie + w^ 1 . Substitution of these expressions into eq. (13. 13j) leads to formula (15. 9p . 

These formulas demonstrate that for Wd, w p 3> w r the effect of measurement shows itself 
in the inverse Zeno effect, in which the effect increases as w r is increased. 

2. Anomalous distribution of measurements 

As in the case of dynamic systems, in our discussion of the anomalous RDM effect in 
the presence of relaxation we will also restrict ourselves to the analysis of the limit of short 
characteristic time between measurements w~ l — > 0, corresponding to ||f2 e ||/uv 1. 
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Most clearly the effect of relaxation can be demonstrated in the limit of fast dephasing, 
when w p 3> v. In this limit the general kinetic equation (15. ip reduces to the system of 
balance equations, i.e. equations for populations of states. 

For the considered two level system (14.11) with relaxation superoperator (15 .3p the balance 
equations for the vector p(£) of populations of states |1) and |2), can be written as 

t> = -R d p, where p(t) = Pl (t)|ll) +p 2 (t)|22), (5.10) 

and 

Ai = tSS(|ll>-|22»((ll|-(22|) (5.11) 

with 

w° d = w d + 2v 2 w p /(w 2 p + 4e 2 ). (5.12) 

The function under study p(e) = Poo{t), which is the Laplace transform of the probability 
Poo(i) = pi(t) of survival in the state |1), can be calculated with the use of expression (14. 161) 
but with supermatrix Q e replaced by Cl e defined in the reduced space of diagonal elements 
of the density matrix 

Poo(e) w (lll^ll^/ail^lll) (5.13) 

. ± ( e + 2 ^) a-1 K 14^ 

e" + (e + 2w° d )« ' l ° ' 

where Cl d = e + R d - 

This formula predicts the long tailed behavior Poo(t) ~ l/t a similar to that found in Sec. 
IVD for dynamic systems. In the limit a — > 1, however, the amplitude of the tail becomes 
negligibly small and the kinetics reduces to the evident exponential: Poait) = e _u, °'. 

Similar to the dynamic systems, the most important feature of the anomalous kinetics 
Poo(t) in the presence of relaxation consists in the independence of this kinetics of w r , i.e. 
of the characteristic time w~ x between measurements. It is worth noting, though, that this 
independence results from the existence of the nontrivial special limit w r — > oo, manifesting 
itself in eq. (I4.16p . rather than from the effect of relaxation itself. 

VI. DISCUSSION AND CONCLUSIONS 

In this work we have discussed the RDM effect on the evolution of quantum systems. The 
sequence of measurements is described as a renewal stochastic process 0, Q] , whose specific 
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properties are controlled by the PDF W(t) of time intervals t between measurements. The 
specific features of the RDM effect, which can be called the stochastic Zeno effect, are 
essentially determined by the analytical behavior of W(t) and properties of the quantum 
system under study, for example the presence of relaxation or decay (see below). 

The analysis carried out in our work shows that in the case of rapidly decreasing PDF 
W(t) the RDM effect on dynamic systems (without relaxation) is similar to that of mea- 
surements equidistantly distributed in time |l|-|4j . This effect is characterized by the average 
time i = J£°dttW(t) between measurements or the average rate w r ~ 1/i of the repetition 
of measurements. To demonstrate the specific features of the effect of the RDM with rapidly 
decreasing W(t) on dynamic systems, we have analyzed the effect on the two level system 
within the Poissonian model for the RDM. The model is shown to allow one to describe both 
quantum and inverse Zeno regimes of the RDM effect by very simple analytical expressions. 

In particular, it is possible to easily treat the main manifestation of the quantum Zeno 
effect consisting in the reduction of the decay rate of the measured state with the decrease 
of the time between measurements. It is also shown that in the limit of frequently repeated 
measurements w r ^> v the RDM result in the exponential decay of the survival probability: 
p{t) ~ e~ Wot [eq. ( 14.13)) ]. with the rate Wo(w r ) which is either decreasing or increasing 
function of w r depending on the relation between w r and the splitting of levels e. These 
two type of behavior correspond to above-mentioned quantum and inverse Zeno regimes, 
respectively 3I. 4} llQ] . 

Of special interest is the anomalous case of heavy tailed PDF W(t) ~ 1 / (w r t) 1+a , (a < 1), 
in which the RDM effect appears to strongly differ from that for rapidly decreasing PDF. 
Note that in this case the average time t does not exist, but the PDF is still characterized 
by the specific time t r = w~ x . For anomalous PDFs the Zeno effect is not observed: in 
the limit t r — > the survival probability p(t) = Poo(t) turns out to be a function of time t. 
Moreover the probability Poo(t) is a very slowly decreasing function: Poo(t) ~ l/t a . 

The relaxation in the quantum system can strongly show itself in the effect of the RDM 
on the evolution of the system. In particular, in the case of rapidly decreasing W(t), for not 
very large rates w r smaller than the characteristic relaxation rate w° d (see Sec. V B), the 
survival probability pit) is demonstrated to exponentially decrease pit) ~ e - """ 4 [eq. (15. 8p ] 
with the rate w®{w r ) either decreasing or increasing as w r is increased. The form of w^(w r ) 
behavior depends on the value of (w r + w p )/e, where w p is the dephasing rate. Similar to the 
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pure dynamic case (in the absence of relaxation) discussed above these types of dependences 
can be considered as a manifestation of quantum and inverse Zeno effects SflQ. The 
acceleration of p(t) decay by measurements, i.e. the inverse Zeno effect, is also found in the 
limit of very fast relaxation u> d , w p 3> w r [eq. (15. 9p ]. 

In the case of anomalous W(t) the relaxation in the system does not lead to any new 
specific features of the evolution kinetics of the system under study. Similar to the dynamic 
systems, for systems with relaxation in the limit of high characteristic rate w r — > oo the 
survival probability Poo(t) is still a function of time slowly decreasing at large times: Poo(i) ~ 



In the end of this section we would like to discuss some recent papers [2CH22| concerning 
the theoretical analysis of quantum measurements and, in particular, the quantum Zeno 
effect within the approaches which has something in common with those applied in our 
work. In above-mentioned papers the kinetics of measurement induced jump-like relaxation 
transitions in quantum systems under study are described both numeric ally by direct Monte 
Carlo modeling of stochastic quantum relaxation transition processes 20] and analytically 
by solving some Markovian kinetic equations 2jJ, l22| for corresponding PDFs. As far as 
analytical approaches are concerned of particular interest is the analysis of Poissonian-like 
iheoretical models of stochastic measurement induced jumps between states made in ref. 



2l| | . It is worth noting though that the authors of this work restricted themselves only to 



the most general (and at some respects too formal) analysis of specific features of Poissonian 
jump processes, including the equivalent description of the process in terms of the SLE or 
corresponding Ito-type stochastic differential equations. Interesting extension of the Marko- 
vian kinetic description of the problem within a quantum model of the measuring device is 



considered in ref 



22(. The perturbative treatment of the system-device interaction allowed 



the authors to derive relatively simple Markovian (Bloch-type) kinetic equations which de- 
scribe the evolution of the measurement affected system. Results obtained in this work are, 
to some extent, related to some results of our study, concerning the Poissonian RDM (since, 
as we have mentioned above, the Markovian processes are directed related to Poissonian 
ones), but unfortunately the direct comparison is hardly possible because of difference in 
model parameters used in both treatments. 

The important contribution of our work, as compared to those discussed above, consists 
in the proposed new approaches, which enable us to significantly generalize the analysis of 
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measurement effect on evolution of quantum systems, allowing for the consideration of the 
RDM effects for different types of stochastic RDM processes (described within the RA) from 
conventional Poissonian one to anomalous. In addition, even in the Poissonian approach, 
which is related to some Markovian models applied in early works, the analysis is essentially 
extended by considering more thoroughly the dephasing effect of measurements and by using 
more general relaxation model for investigating the manifestation of relaxation in the system. 

Concluding our short analysis of obtained results we would like to discuss the possibilities 
of experimental observation of the effects predicted in our work. Of course, the question of 
great importance is whether it is possible to realize the RDM experimentally. Answering 
this question it is worth emphasizing the following points: 

1 . The RDM can be realized experimentally with the use of equipment which can random- 
ize the set of measurements. Among possible variants of random processes the Poissonian 
one is, probably, the most simple for realization. At first sight the corresponding efforts will 
be of not very much use in reality, however, it is worth noting that thus obtained results on 
RDM effect can be analyzed somewhat easier than those of equidistantly distributed mea- 
surements. The simplicity results from the possibility of the analysis with simple Markovian 
kinetic equations for any number of experiments. In addition, as it was mentioned in the 
above discussion, the Poissonian RDM can be realized in the measurement processes them- 
selves due to stochastic nature of the measurement procedure in some particular experiments 
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2. Perhaps, the most important fact is that the RDM can be realized by the process under 
study itself. The point is that in many cases the measurement is made within fairly small 
volume which can be a small part of the larger volume where the process occurs. For example, 
suppose that the measured system is a small Brownian particle with quantum internal 
degrees of freedom, undergoing stochastic motion within large volume. In its motion the 
particles crosses the small measurement volume. Any visit of this volume can be assumed to 
result in the measurement. In such a scenario of measurements their statistics is represented 
by that of visits of the small measurement site. 

The statistics of visits depends on the mechanism of motion of the particle. For instance, 
in the case of the particle, confined in relatively small volume (of type of a cage) the statistics 
is close to Poissonian. For freely diffusing motion in the infinite space the statistics is quite 
well described by the anomalous variant of the RA with the heavy tailed PDF W(t) ~ l/t 1+a , 
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in which a < 1 is determined by the dimensionality d s of the space js, 9]. 

3. As a possible modification of the procedure, described above in p. 2, one can consider 
the "measurement" of the quantum state of a Brownian particle by another large quencher 
particle whose effect on the quantum subsystem of the migrating Brownian particle can be 
treated as a measuring device. The statistics of measurements in such systems is, actually, 
determined that of reencounters of the Brownian particle with the measuring quencher. The 
properties of this statistics is well described by the RA, as mentioned above, with the specific 
features depending on the mechanism of migration. 
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Figure captions. 

Fig. 1. The dependence of the survival probability p(t) = p(r|r r ) on the inverse average 
time between measurements r r = w' 1 = v/w r for two dimensionless times r = tv [v is 
defined in eq. (14.21) ]: r = 5 (a) and r = 10 (b). The probability is calculated for the two 
level system (14. 2p . (14.3)) . assuming the Poissonian statistics of measurements. Calculation is 
made with the use of formula (14. lip for four values of e = e/v: e = 2.5 (triangles), e = 1.0 
(circles), e = 0.5 (squares), e = 0.25 (stars). 

Fig. 2. The dependence of the survival probability Poo(t) on r = tv, calculated by the 
inverse Laplace transformation of Poo( e ) (14.16)) . for anomalous RDM with different parame- 
ters e = e/v and a: (a) e = 0.94, a = 0.1 (full line); e = 0.53, a = 0.3 (dashed line); and (b) 
e = 0.53, a = 0.92 (1), e = 0.94, a = 0.92 (2); e = 0.53, a = 0.97 (3); £ = 0.94, a = 0.97 
(4). In Fig. 2a the straight lines corresponding to the dependence Poo(t) = 1/(2. 3t a ) are 
presented for the sake of demonstration of the asymptotic behavior of the exact dependence 
Pooit). In fig. 2b the straight (dashed) lines correspond to the approximate dependence 
(Q9|) . 
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